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We consider the disordered quantum rotor model in the presence of a magnetic field. We analyze 
the transport properties in the vicinity of the multicritical point between the superconductor, phase 
glass and paramagnetic phases. We find that the magnetic field leaves metallic transport of bosons 
in the glassy phase in tact. In the vicinity of the vicinity of the superconductivity-to-Bose metal 
transition, the resistitivy turns on as (H — H c ) 2 with H c . This functional form is in excellent 
agreement with the experimentally observed turn-on of the resistivity in the metallic state in MoGe, 
namely R w R C (H — -Hc) M , 1 < fJ, < 3. The metallic state is also shown to presist in three spatial 
dimensions. In addition, we also show that the metallic state remains intact in the presence of Ohmic 
dissipation in spite of recent claims to the contrary. As the phase glass in d — 3 is identical to the 
vortex glass, we conclude that the vortex glass is, in actuality, a metal rather than a superconductor 
at T — 0. Our analysis unifies the recent experiments on vortex glass systems in which the linear 
resistivity remained non-zero below the putative vortex glass transition and the experiments on 
thin films in which a metallic phase has been observed to disrupt the direct transition from a 
superconductor to an insulator. 



I. INTRODUCTION 

When superconductivity is quenched in thin metal alloy films, a metallic 1,2,3,4 rather than an insulating phase 
ensues. The subsequent metallic phase is robust as it exists over a wide range of magnetic field and film thickness 
(disorder). Regardless of the statistics of the charge carriers, a T = metallic phase in two dimensions does not mesh 
easily with conventional wisdom. In the case of bosons, the uncertainty principle dictates that in the absence of some 
new state, bosons either condense or localize as in conventional treatments 5 of the insulator-superconductor problem. 
Further, non-interacting electrons^ carry no curent in two dimensions in the presence of disorder. Precisely what 
the new exotic metallic state is is hotly debated 7 ^. While many proposals have been made 7 -, only the phase glass 
model 7, has been demonstrated through direct calculation to have a finite resistivity at T = in the presence 
of disorder. In this model, the charge carriers are strictly bosonic. Metallic transport arises from the coupling of 
the bosonic degrees of freedom to the low-lying excitations of the glassy phase. The densitj*2*ii*i2*i 3 " of such low-lying 
excitations scales as \oj\ and hence gives rise to a dynamical exponent of z — 2. Indeed, other models have been 
proposed which could give rise to z = 2 dynamics. For example, a popular view is that the metallic state for the 
bosons is driven by dissipation arising from normal electronic excitations^ 4 * 1 ^. Such fermionic excitations are typically 
modeled by a dissipative Ohmic heat bath. However, bosons described by an Ohmic propagator of the form, 

G Q = (k 2 + V \LJ n \+m 2 )-\ (1) 

still only exhibit either superconducting or insulating phases^i 7 -. Aside from providing a dissipative environment, 
normal electronic excitations could carry the current in the metallic state. As the normal state of the materials of 
interest, MoGe for example, is well described by Fermi liquid theory, the localization principle, which precludes a 
metallic state for non-interacting electrons in two dimensions^, should be operative. Hence, while normal electronic 
excitations are undoubtedly present above the superconducting state, invoking— them to explain the origin of the 
metallic state is a non-starter. Additional fermionic models have been suggested^ in which bosonic vortices are 
attached to a strongly fluctuating gauge field which transmutates the statistics. On phenomenological grounds, the 
resultant fermionic vortices were assumed to diffuse and hence constitute a gapless liquid. While flux attachment 
schemes have been highly successful in describing gapped phases 2 ^, such threadings are problematic in gapless states 21 - 
as the fluctuations of the gauge field are relevant. 

In the phase glass model, fermionic excitations need not be invoked because the dissipation is self-generated from 
the glassy environment. The resultant bosonic propagator (derived in the Appendix A) is of the form, 

Gf b \k,uj n ) = G (k, uj n )S ab + /9Gg(k,a; w )«J Wwi o (2) 

where a and b are the replica indices, q is the Edwards- Anderson order parameter for the glassy phase and u„ are the 
Matsubara frequencies. The first term is the standard bosonic propagator. However, the second term arises entirely 
from the glassy landscape. Although Eq. J5J) can be derived explicitly, it also follows directly from dimensional 
reduction. That the second term contains the square of the bare bosonic propagator, Zinn- Justin^ surmised holds 
profound consequences for disordered systems. Indeed this is true. It is from this term that the metallic state arises. 

To constitute a complete theory of the metallic state that disrupts the insulator-superconductor transition, the 
phase glass must be applicable to systems in the presence of a magnetic field as well. It is precisely this question that 
we address here. We show that the metallic state persists in the presence of a magnetic field. In the weak-field limit, 
we find that the dominant role of the field is to renormalize the inverse correlation length. The low-field limit of the 
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resistivity is described by a power law of the form R C {H — H c ) 2zL/ , where z = 2 and v = 1/2 at the mean-field level. 
This functional form is in excellent agreement with the experimentally observed turn-on of the resistivity in MoGe 2 ^. 
We also find that the metallic state survives in 3D. In 3D the relevant comparison is with the vortex glass^ which 
has been argued extensively to be a superconducting phase. We find instead that this state is metallio^iiSSiSi. 



1. Phase Glass Overview 



We begin by summarizing the key elements of the zero-field theory we will need when we consider the effects of a 
non-zero magnetic field. The phase glass arises naturally in the context of the disordered quantum rotor model 



(i,3) 



with random Joscphson couplings Jy but fixed on-site energies, Eq- The phase of each superconducting island is 
6i. Note that additional on-site disorder of the form ivjd/dOj results in the equivalent particle- hole symmetric field 
theory provided that the distribution of on-site energies has zero mean. For the on-site energies, the non-zero mean 
case is irrelevant here as this corresponds to a density-driven insulator-superconductor transition (IST)2&. 

If the Joscphson couplings are chosen from a distribution with zero mean, only two phases are possible: 1) a glass 
arising from the distribtion of positive and negative J-jS and 2) a disordered paramagnetic state. A superconducting 
phase obtains if the distribution 



1 

V2^J2 



P(Jij) = 7= Tn exp 



2J 5 



(4) 



of J-jS has non-zero mean, Jo, and J the variance. To distinguish between the phases, it is expedient to introduce the 
set of variables Sj = (cos 0i, sin 9i) which allows us to recast the interaction term in the random Josephson Hamiltonian 
as a spin problem with random magnetic interactions, Y^u j\ JijSi ■ Sj. Let (...) and [...] represent averages over the 
thermal degrees of freedom and over the disorder, respectively. Integrating over the random interactions will introduce 
two auxiliary fields 

Q^(k,k',r,r')-(^(k,r)5^(k',r')) (5) 

and $"(k,r) = (S?(k, t)), respectively. The superscripts represent the replica indices. A non-zero value of ^"(k, r) 
implies phase ordering of the charge 2e degrees of freedom. For quantum spin glasses, it is the diagonal elements 
of the Q-matrix D(t — r') = lim Jw o jg^(Q^(k, k', r, r')) in the limit that |r — r'| — * oo that serve as the effective 
Edwards- Anderson spin-glass order paramete»ii*i2ii2i within Landau theory. The essential aspect of the quantum 
rotor spin glass is that the saddle point solution for the corresponding action is minimied by a solution of the form 

Q#(k,wi,wa) - p{2Ti) d 5 d {k)8^[D^i)^ 1+lJj2fl 5 ab 

+/3<Wo<W<?° 6 ] (6) 

where the diagonal elements are given by 



D(u) = -Voj 2 + A 2 /k, (7) 

with k a coupling constant in the Landau free energy for the spin glass. The diagonal elements of the Q-matrices 
describe the excitation spectrum. Throughout the glassy phase, A = and hence the spectrum is ungapped and given 
by D(lu) = —\lu\/k. The linear dependence on arises because the correlation function Q„u(r) decay o 11 ' 12 ^? as t~ 2 . 
This dependence results in a fundamental change in the dynamical critical exponent from z = 1 to z = 2. In a recent 
treatment of the disordered quantum rotor problem, Ikeda™ included an Ohmic dissipative \u>\ term at the level of 
the Hamiltonian. He has found that the excitation spectrum in the spin-glass order parameter scales as M 1 / 2 . With 
the M 1 / 2 ansatz, Ikeda 2 ^ has found that the resultant phase is a superconductor. We shown in Appendix B, however, 
that this claim is without merit. Namely, even if \u>\ dissipation is put in by hand at the outset, the universality class 
of the problem does not change. Such a contribution simply renormalizes the coefficient of |cj|. 
To analyze the bosonic conductivity, we focus^ on the part of the free energy 

rm= (fc 2 +^+™ 2 )i^(k,^)i 2 

a,fi,k,uf n 
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-1 J d d X J d Tl dT 2 ^(x,T 1 )^ b u (x,T 2 )Q a Jl(x,T 1 ,T 2 ) 

+U [ dTj2[^(x,T)*°(x,T)] 2 (8) 

governed by the fluctuations of the superconducting order parameter. The coupling between the spin-glass and bosonic 
degrees of freedom fundamentally changes the dynamics as can be seen by substituting the saddle-point solution for 
the Q- matrices into Eq. JSJ). The resulting expression for the free energy 

•fgauss = (k 2 +U 2 n + V \LU n \+m 2 )\lp a {k,LU n )\ 2 

-f3q ^ ^oVAk^^k,^)]*, (9) 

reflects the damped dynamics introduced by the glass degrees of freedom. As remarked earlier, the |w| term dominates 
giving rise to z — 2 dynamics. Certainly, dissipative models with \uj\ dynamics^ put in by hand have been studied 
explicitly in the context of the insulator-superconductor transition partly in the hope of finding new phases 1 - 5 . However, 
direct calculations^ of the conductivity reveal that \oj\ dynamics alone cannot give rise to a metallic phase. What is 
new here is that the \lu\ term is not put in by hand but arises naturally from the coupling between the glassy dynamics 
of the phase and the charge degrees of freedom. Further, as is evident from Eq. J5J), the Gaussian propagator changes 
fundamentally as a result of the ulrametric properties of the glassy fields, g ab . Substitution of Eq. (0 into the Kubo 
formula for the conduct ivityS 

^ = 0,2^0) = ^ (10) 

demonstrates that both the \uj\ dynamics and the glassy landscape conspire to yield a non-zero bosonic conductivity 
at T = 0. This expression crosses over smoothly to the a — oo conductivity in the superconducting phase (to = 0) to 
the vanishing conductivity of the insulating state (q — 0). Hence, the quantum rotor model in which the disordered 
Josephson couplings possess a non-zero mean describes a metal, insulator and a superconductor, 



2. Magnetic Field: Summary of Findings 



In this paper, we extend the formalism summarised above to include a magnetic field. Since disorder is also included, 
our results in 3D are relevant to the vortex glass stated. In the weak field limit, we find that the transport properties 
are well approximated by replacing to 2 in the conductivity by 

to 2 = to 2 + ^m 2 H (11) 

where m 2 H = 2e h H — 2/£ 2 is the inverse square of the magnetic length (£) in units of c = 1. Physically, this result 
is transparent as the predominant role of the magnetic field is to provide a confining potential for the electrons on a 
scale set by 1/yJ ra 2 H . The new condition for the onset of superconductivity is 

m 2 + 2 TO ff = superconductivity criterion (12) 

which has a non-trivial solution because to 2 < demarcates the superconducting phase. Hence, we find quite generally 
that metallic transport survives in the presence of a magnetic field. Further, we find that metallic transport persists 
even in d=3 where the vortex glass phase has been proposed to exist in the presence of disorder—. Because of 
the equivalence of the phase glass to a vortex glass in the presence of magnetic field, our results imply that the 
linear conductivity should remain finite in the vortex glass. These results are discussed extensively in the context of 
the insulator-superconductor transition and the recent experiments that indicate that the vortex glass state is not 
accompanied 2 ^**^ by a vanishing of the linear I-V characteristics. 



II. PHASE GLASS CONDUCTIVITY IN A MAGNETIC FIELD 



In the presence of a magnetic field, the 0(2) quantum rotor model becomes 

' d_ 

86. 
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where Ay = (e*/h) J? A • c?l (e* = 2e). The Josephson couplings are assumed to be random and governed by the 
distribution Eq. In the current work, we treat A to be a static time-independent magnetic field as this is the 
most experimentally relevant situation^. For a random system, the technique for treating disorder is now standard^: 
1) replicate the partition function, 2) perform the average over disorder and 3) introduce the appropriate fields to 
decouple the interacting terms that arise. As the corresponding action has been detailed previouslyii^, we will 
provide additional steps that are necessary to determine how the electromagnetic gauge couples to the spin glass 
order parameter. We write the replicated partition function as 

v.jven,^ (i4) 

where the Euclidean action is given by 



s = E^(1r) 2 -EE^»K a w-^M-^w] ( 15 ) 

J0 { i,a ° V 7 a (ij) ) 

where the superscript a represents the replica index. The integration over in Eq. I|14(l results in the effective 
action, 



+ 1tEE/ f dTdT '\ £ e X p{z[^(T)-a^(T')-(^(r)-a^(r'))-(^(T)-a^(rO)]} 

a,b (ij) J ° J ° a=+l.-l 

+ C.C. 

-yEE f dre^{9t-e){r)-A l3 )+c.c. (16) 

a ij 



with a = +1, — 1. As a result of the sum over a, we see that the vector potential enters both symmetrically and anti- 
symmetrically. As in the zero- field case, we introduce the two-component spinor, S a (r) and the auxiliary field defined 
in Eq. © to describe the spin-glass transition. The remaining steps involve performing the cumulant expansion and 
taking the continuum limit. The final action can be separated into the local and gradient parts: 

S e « = Sloe + S gI (17) 

where the local part 



>loc 



d d 



-~ J dr 1 dT 2 dr3^Q^(x,r 1 ,T 2 )Qtp(x,T 2 ,r 3 )g-(x,r 3 ,r 1 ) 

+ J d d x l - J rfr^[ U Q^(x,r,r)g-( X ,T,T)+«Q^(x,T,r)Q^( X ,T,r)]| 
-§ / d d x J dT.dT^iQ^r^T,)] 4 

a , b 

+ Jd d xJ dr^-^^(x ) T 1 )-^^(x,T 2 )\ Tx ^ 2 = T + ^ J dr t^^t) 

a " a 

--^J d d x J dndr 2 X)^(^n)^(^r2)Q^(a ; ,T 1 ,T2) + --- 



(18) 
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contains a spin-glass part that is identical to that of Read, Sachdev and Yeii , a superconducting part derived 
previously and a gradient part 



d d x / dr\dr2 



if if 
V - — A(x, n) + — A(x, r 2 ) ) Q a _ b (x, n, r 2 ; 



d d 2 



d,Tld,T2 



a . 6 



V- ^-A(x,ri) - ^-A(x,t 2 ) ) r/'"(x. -,. 



-7- ( V-— 



(19) 



in which the vector potential couples both symmetrically and asymmetrically to combinations of the Q— matrices 
of the same parity as well as to the bosonic part of the order parameter. Using the fact that Q±(ti, t 2 ) ~ 
(exp [i(0f(r] ± ^(t'))]), we write the parity combinations of the Q— matrices as follows 



QS?(x, n.Ta) = \ [Qf 1 (K,r 1 ,T 2 )TQ^,ri,r 2 )]+^[Qf 2 (^n,T2)±Qf 1 (^r 1 ,r 2 )] 



(20) 



It is evident that the vector potential enters in a non-time translationally invariant manner. This is a direct conse- 
quence of the fact that the Q— matrices themselves are a function of two independent times, not simply the difference 
of ti - t 2 . 



1. Transport in 2D 

Our immediate focus is on the bosonic fluctuation conductivity in the glass phase. To this end, it suffices to focus 
on the i/j-dependent part of the action: 



V-^U(*,T 

en 



*2(*,r) 



+ jd d xj dr^^l(x,n)^l(x,r 2 )\ Tl =r 2 =r + ^ f drJ2K(^r)%( Xl r)] 2 

J a 1 2 a 

I d<lx I d nrfT 2 ^^(^ ! r 1 )^(a ; ,r 2 )Q^(a ; ,T 1 ,T 2 ). (21) 

a, b 

A precursor to the final result can be obtained by assuming that ij) is independent of both space and time. In this case, 
the only vector-potential dependent terms that survive in the total action are those arising from the pure diamagnetic 
term, A 2 , in the bosonic part of the action. For the field along the z— axis, we write the vector potential as 



A=-(Bx r) 

2 V ' 



(22) 



Consequently, the integral of the diamagnetic term over a rectangle with dimensions L\ x L 2 yields (L 2 + L 2 )/3£ 2 . 
This result can be combined into the gradient part of the action for the bosonic degrees of freedom simply by redifining 
the mass term as 



l^l+-(L 2 + L 2 )/£ 2 = ^. 



(23) 



At 



field 



V> 2 = ^(7-A) 



(24) 



provides a non-trivial phase boundary for the superconducting state. In the presence of a magnetic field, our prelim- 
inary considerations yield to a modified boundary 



7 > A. 



(25) 
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Because 7 > 7, we find that the presence of magnetic field suppresses the superconducting region. Because 7 acts as 
the mass term for the superconducting order, an intimate connection 



2 Jnd - E, 



c 



2J da 2 



(26) 



exists with the basic coupling constants, Ec and Jo- Here we have adopted the mean-field expression^ with a is the 
distance between the superconducting islands and d is the spatial dimension. Use of 7 in this expression results in 
a mean-field expression for the critical field. Hence, the renormalization of the mass by the inverse magnetic length 
has immediate physical consequences. For example, by simply replacing m by 7, we obtain the leading magnetic-field 
dependence of the conductivity. A more rigorous treatment will show that this result is almost correct. 

A correct treatment of the conductivity requires that we relax the constraint that tp be spatially independent. In 
this case, we need to expand ip ina basis of Landau levels. To this end, we define 



— pVyV 



Mx)C a (l, P y) 



(27) 



where the (j)i( x ) are the Hermite polynomials. Within the Landau gauge, the gradient part of the action requires the 
solution to a Landau level problem of the form, 



dx 2 



dy 2 



1 



d 



2i 2Py dx 



1 



1 + 77 m H 4>i(x). 



(28) 



Expanding the free energy in the Landau level basis results in the replacement of momentum summations with sums 
over Landau levels. The resultant free energy 



Fib 



E 

Um 2 H 



(I 



1 , 



'H 



47T/3 



dr 



E 



//|^„,| \C a (l,Py, L0 m )\ 2 

\C(l,p y ,u! m )\ 4 + l3q8 UJm fiC a *(n,py,uj m )C h (l,Py,uj m ) 



(29) 



can be written entirely in terms of the expansion coefficients C a (l,p y ,Lu m ) once the Landau level problem is solved. 
Our expression for the free energy lays plain that the naive association of the mass with m 2 + m 2 H /2 is not entirely 
correct as an Z-dependent term exists as well. The corresponding Gaussian propagator is of the form, 



Gtf(l,u> n ) = G {l,w n )6ab + /3G^l,u n )q5 Ua ,o 



(30) 



in the n — > (zero replica) limit with Go(l, tu n ) = (m^(i + 1/2) + uj 2 + t]\ui\ + to 2 ) -1 . To calculate the conductivity, 
we add to the vector potential a small perturbation 5A(r, t) which we use as a source to calculate the current and the 
subsequent conductivity which according to the Kubo formula 



a(iu n ) 



lim — 



d d (x-x') / d(r- 







S 2 Z- 



'SA x (k, t)6A x (x>,t>) 



(31) 



involves two derivatives of the free energy with respect to the vector potential. The resultant Kubo formula for the 
conductivity^^ 



2h@u) n 



j2 j2^ i+i n 2G t( i ^ m )G a b (i+i,o Jm ) 
=0 



-Gg b (Z, tu m )Gf(l + 1, uo m + oj n ) - Gf(l, cj m + cj n )G a b {l + 1, w m )] 



(32) 



involves an explicit sum, I, over the Landau levels. In the absence of the glassy landscape (that is, q ab = 0), this 
expression has been evaluated previously 35 36 . The result in the quantum disordered phas o 35 i 36 is that an insulating 
phase obtains regardless of the magnitde of the dissipative dynamics \u>\. Consequently, the only terms that might 
give rise to a metallic phase are those proporotional to the Edwards- Anderson order parameter, q. The leading such 
term gives rise to a conductivity of the form, 



cj 2D {lj = 0,T ^0) = 



(e*m 2 H ) 2 
h 

(e*m 2 H ) 2 



1=0 



(i + 1) 



2q V 



(m 2 + '±)m 2 H ) 2 {m 2 + (I + \)m 2 H ) 2 



(33) 
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Where x is defined as x = ^^Sr- The summation over I give rise a digamma function defined as ^(1,0;) = dln fJ' x ' > 
where T(x) is the Gamma function. The large and small field limits relative to to 2 can be found exactly. In the weak 
field regime, we find that the dc conductivity 



a 2D [co 



(w = 0,T-fO) 



4e 2 qr\ 



<H 



4e 2 qrj 



7 

20' 



x — mfj/rn 2 <C 1 



(34) 



can be written compactly in terms of the renormalized mass scale, to. In obtaining this expression, we worked in the 
correct order of limits^! for the dc conductivity; that is, u> = 0, T — > 0. When the field vanishes, this conductivity 
goes over smoothly to Eq. (|10|) . In the opposite limit, m 2 H ^> m 2 , the conductivity reduces to a simple sum 



<7 2D (uj = 0,T 



0) = 



(e*) 2 qq 



•h 



2 + — 

x 



X > 1 



(35) 



Consequently, we find then that the conductivty of the Bose metal remains finite in both the weak and strong field 
regimes. 



2. Transport in 3D 

That the fluctuation conductivity reaches a plateau in the phase glass in the presence of a magnetic field does not 
appear to have been expected previously. Before we consider the effect of the quartic term, we show that the metallic 
state persists in three spatial dimensions as well. The key in generalizing the formalism presented thus far to 3D is 
to include in the expansion for the order parameter, 

r(l,x,y,z, Py ,p z ) = e^y+^ z Mx)C a (l,Py, Pz ) (36) 

a plane wave for the third direction. For a vector potential of the form, A(r) = (0, Hx, 0), the spatially varying part 
of the action involves the eigenvalue problem, 



d 2 , , , . fe*H 
dx 



2 ipi{x)+[—-x Q ) i>t(x) =E l {p z )^ l {x) (37) 



where 



E = m 2 H {l+ l -)+p 2 z . (38) 

Except for the p z -dependence, this oscillator problem is identical to the 2D analogue. In fact, expanding the free 
energy in terms of the states in Eq. I|3t)|) and computing the derivatives with respect to the vector potential in the 
Kubo formula reveals that the primary difference in the longitudinal conductivity 

a* D {iu> n ) = [S ]T ^(/ + l)[2Gg b (/,c m ,p z )Gg b G + l, Wm ,p 2 ) 

a,b,pz,0Jm t— u 

-Gf(l,uj m ,p z )Gt(l + l,w m + u n)Pz ) - G a b (l,Lu m +uj n ,p z )G a b (l + l,u m ,p»)] (39) 
is the presence of a sum over p z . Because the p z dependence in the bare Green functrion 

Go(l,p z ,w n ) = (lm 2 H +ul + i 1 \uj\+m 2 +p 2 z y 1 (40) 
amounts to a shift in to 2 , the final expression for the conductivity in D = 3 

^(^0,T^0) J e 'fV f(l + l) n 2 „ 0/0 , 2 2 ^ 2 2 -- 2 , 2 . (41) 

h ( m + +Pz + ( l + 3/2)to^) 2 (to 2 +p 2 z + (I + l/2)m 2 H ) 2 

is easily obtained from Eq. I)33|l . Let us define 

2 ( 1' 



An = m z + ( n+ - ) m- H . (42) 
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FIG. 1: Conductivity of in the presence of magnetic field. In 2D, the longitudinal conductivity is given by a 2D = 0x0/2(3;) + 
o"iff2(a;) with 00 = 8e 2 qrj/hm 4 and a\ = 8e 2 q 2 Ur]/hiTm (i . The corresponding expression in 3D is a 3D — maofs(x) + mcrigzix). 



Performing the integral over p z using residues reduces the conductivity 



a(uj = Q,T -> 0) = 



2% 



4X>+i) 



71=0 



1 



Aj Aj +l {Aj + A^ +1 ) A n A n+ i(A n ' + Aj +1 ) 3 _ 



,1/2 



1/2 



,1/2 



(43) 



to a single sum over the Landau index. As in the 2D case, this expression can be evaluated analytically only in the 
small and large field limits. Defining x = m 2 H /m 2 1 we find that the strong-field limit is given by 



a 3D = ( e *) 2 QW_ 

h to 3 

where we have made the approximation, 



x — m 2 H /m 2 1, 



(44) 



E 



(n + l/2)™(n + 3/2) 1 /2 



2m + 1 



(45) 



In 3D, the leading asymptotic inverse correlation length dependence scales as to 3 as opposed to the m dependence 
obtained in 2D. In the opposite limit, 



r 3D 



(e*) 2 qyir 

h m 3 



1 1 

- + -x 

6 8 



1. 



(46) 



For the sake of comparison, we have evaluated the conductivities in 2D and 3D numerically and plotted them in Fig. 
as a function of x. As is observed experimentally, both are decreasing functions of the magnetic field. 



Role of interactions: 



Thus far, we have worked at the Gaussian level in which interactions are strictly ignored. At the tree level, a 
dynamical exponent of z = 2 renders the quartic interaction U marginally irrelevant. However, considering the last 
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n-l,Q, 




n+1, co m 

FIG. 2: One-loop diagrams that determine the leading interaction corrections to the conductivity. The squiggly line represents 
the quartic intera ction . The matsubara frequencies and Landau level indices are indicated. These determine the conductivity 
correction in Eq. 1471 . 

term in Eq. JJjJ on equal footing with U in the one-loop renormalization group scheme, we reach the conclusion that 
the RG equations flow to strong coupling. The relevance of q at all dimensions manifests itself also by the increasing 
singularity of relevant contributions from higher order diagrams in the perturbation series in U . Consequently, we 
consider the role of interactions. Consider first the linear U correction which is given by the diagrams in Fig. J5J. 
We use the definition of A n to simplify the notation in the evaluation of the diagrams. In two dimensions, the sum 
of these diagrams is given by 



(1) _ (e*) 2 q 2 Um 2 H 



2fo p u ^2^^(n+ l)(j4 n ) 7i A n ] rl + A n ^ 1 A n 1 + 2A n 2 A n 2 1 ) 
- £)(n + l)m% (2A-l 1 G 2 (n, u n ) + 4A- 3 G (n + 1, u n )) 

n 

-5> + lK, (2A- 2 G 2 Q (n + l^ n )+4A-l 1 G Q (n,u Jn )) 1 £ Af, 

n J I 



which simplifies in the zero-frequency limit to 

(i) _ 32e 2 q 2 Ur] ^ (n+l)x 3 



E 



h 4?rm 6 A 3 n A 2 n+iA 2 ■ 



(47) 



(48) 



Two relevant features are that the leading U correction has a non-divergent contribution to the conductivity in the 
T — ► limit and the dependence on m is more singular than in the U = case. This points to a scaling function of the 
form a « (e 2 /%){rjq/m A )^ (q/m 2 ) f(U,u)u) where <£>(?/) ~ y p for large y, which yields the critical behavior a ~ mT x 
with x = 4 + 2p and f(U, ojh) is a general some function of the interactions and the magnetic field. The value of the 
exponent p cannot be inferred at any finite order in perturbation theory. 

The analogous expression can be derived for the three-dimensional case. The only complication here is the integra- 
tion overp z . The result, 



<J3u(un = 0) 



4e 2 q 2 r)U 



E 



3tt n + 1 7T (n + 1 ) (A n A 1 ^ 
8 A 5/2 + 2 

sin 



OA 1 / 2 A 



.3/2.3/2 M l/2 



A 1 ' 2 V 



{A n +LU H ){n + l) 



2 x 3/2 .3/2 M 1/2 .1/2 s A 
sin A n+1 {An + A n+1 ) J A 



1 



,3/2' 



(49) 
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FIG. 3: Disorder contribution to the renormalized gap. In 2D, Sm% D — (qU / '4nm 4 )h,2(x) whereas in 3D, 5rn? iD = 
(qU/4-Km 3 )h 3 (x) . 



expressed in terms of the function A n is once again finite as T — > 0. Hence, the bosonic conductivity is robust to 
interactions (at least in leading order) in both D = 2 and D — 3 in the presence of a magnetic field. 

Of course, the mass, m, is also renormalized by the interactions. To understand its dependence on the interactions 
at large N. The self-consistent 1/N expansion for the correction to the mass gap, 

53 G o(um,n)+Pq5 Um , G (uj m ,n) 2 (50) 



U 4tt 

contains two terms. The term independent of the glass environment has been calculated previously 3 ^. The second 
term, written succinctly as, 

5ml D = uj H Tf3q \ ^ 1 

U 47rm 4 ^ (1 + (n+ l/2)x) 2 { ' 

is explicitly independent of temperature and provides hence only an added magnetic field dependence to the gap 
equation. Here x = m 2 H /m 2 . In the 3D case, the integral over p z leads to slightly different power in the denominator, 



6m lp 



5 (l + (n + l/2V) 3 / 2 ' (52) 



U 8m ^ (l + (n + l/2)a;) 3 / 2 ' 

with, once again, no added temperature dependence. Hence, the interaction terms that couple to the glassy degrees 
of freedom lead to an innocuous renormalization of the mass as depicted in Fig. iJSJ. 



4. Absence of phase stiffness 



In the strict sense, a physical system possesses a non-zero phase rigidity if upon a non-uniform rotation of the 
phase, the free energy increase is of the form, 

AF=$ (<?r\V9\\ (53) 
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where p s is the spin or superfluid stiffness and 9 is the collective phase variable. For the problem at hand, it suffices 
to expand the gradient part of the free energy for the spin-glass part in terms of a local part given by Eq. © and a 
spatially varying part we refer to as Qff v (k,u)i,u)2)- 

Firstly, let us consider the contribution of the spatial independent part of Q. It has been shown recently^! that a 
Q th term which is responsible for replica symmetry breaking in the spin glass action changes the mean-field solution 
for the excitation spectrum, D(u>). In this case, the excitation spectrum, 

D(w) = cq% A - M/k (54) 

remains gapped in the spin glass phase. Here c is a constant and qEA is the diagonal part of q ab . In the absence 
of the Q th term, the stiffness has been shown to vanish^. However, recently, Ikeda^ has disputed this claim and 
has proposed that the constant term in Eq. (|54|l gives leads to a non-zero stiffness. We show that this claim is false. 
Following our earlier calculation^ (which IkedalS also follows), we find that the relevant contribution to the spin-glass 
stiffness containing the D(uS) term reduces to 



16e 2 

<Jps(lLOn) = T 

JVjJn 



Tj2D{cJ m ) (D(uj m ) - D{u m + u n )) - q EA (%) - D(0)) 



(55) 



in the presence of the constant term. As we have shown previously^, ops can be evluated by analytically continuing 
the first term using the definition of D{u>). The result is non-critical and metallic even in the presence of the 
constant contribution to D(uS), This state of affairs arises because as long as difference of D's always appear, any 
constant that is added to D(0) cannot yield a divergent contribution to the spin-glass conductivity, contrary to the 
claims of IkcdaSS. 

The contribution of spatially dependent part of Q can be extracted from the propagator for the spatially varying 
part of the action. Note that from Eq. (J2DJ, Q a } is determined by the longitudinal components of the spin-glass order 
parameter only. Further, Q ab does not couple to the magnetic field. The transvere components enter only and 
couple directly to the field. Hence, it is only the spatial fluctuations of that are of interest here. The associated 
propagator for Q°^(k, u>\, UJ2) at the Gaussian level is simply 

gf{ Jl ,u }x , U2 ) = T — T 1 (56) 

(n + ±)m H + + \uj 2 \ 

Because w is a Matsubara frequency, the propagator describes the motion of mode whose energy dispersion scales as 
i(n + l/2)m 2 H . Such a mode is purely imaginary and hence does not possess a stiffness. For Q"^, (n + l/2)m 2 H is 
replaced by k 2 and diffusive transport ensues. Hence, neither the transverse nor the longitudinal components possess 
a stiffness. Consequently, the phase glass is not a superconductor and only the fluctuation conductivity computed 
previously contributes to the dc transport. A corrolary of this work is that the vortex glass, in which the disorder in 
the Jy's (caused by a random vector potential) is bounded between [—1,1] does not have a stiffness and hence is not 
a superconductor. 



5. Comparison with Experiment 

Our theory makes precise claims as to how the resistivity turns on in both the disordered and magnetic-field 
tuned transitions. Hence, given the plethora of experimental data, it behooves us to put our predictions to the test. 
Consider first the disorder-tuned transition. In this case, the disorder resistivity turns on as the fourth power of the 
mass according to Eq. (fTOft . As m 2 oc and £ cx (d — d c )~ u , we find that the conductivity turns on as (d — d c ) 2zu . 
For the problem at hand, the mean- field exponents are given by z = 2 and v — 1/2. Any attempt to go beyond 
mean- field will result in an uncontrolled expansion as a result of the coupling to the Q ab degrees of freedom— 1 . If a 
controlled expansion were possible, v would be renormalized not z. Hence, to understand how close the transition is 
to mean-field, we treat v as a fitting parameter. Consequently, we used the expression p(T = 0) oc [d — d c ) Au to fit 
the data on three different samples. The results shown in Fig. Q point to a value of v ranging between [1/2,2/3]. 
To corroborate this prediction, we turn to the magnetic field tuned transition. Mason and Kapitulnik 23 find that the 
resistance per square turns on as 

R cx R C (H - H c y 1< (i < 3. (57) 

Since we are interested in the region close to the lower critical field, it suffices to investigate the critical properties 
of the conductivity using the weak- field expression, Eq. I|34|) . A rough scaling of the conductivity can be obtained 
by noting that in the weak-field expression, m 2 is replaced by m 2 + m 2 H = rh 2 . If we replace rh 2 by h~ zu , with 
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log[(d-d c )/d ] 



FIG. 4: Critical scaling of the turn-on of the resistivity near the superconductor-metal transition on the metallic side as a 
functin of the sample thickness. d c represents the critical thickness needed to destroy superconductivity. Based on Eq. (1I0L 
we are able to extract a value of the correlation length exponent, v. The fits establish a value of v roughly between [1/2, 2/3]. 
The data were extracted from Ref. (1). 



h = (H — H c )/H c , we find that p w h Av . To fit Eq. I|57|) . we require that 4i/ ss 2.6 or alternatively, v — .63, a value 
not inconsistent with the disorder-tuned transition, though it is unclear the precise magnitude of the error bars in 
this fitting procedure. Consequently, we resort to a more rigorous treatmet. The transition to the metal state occurs 
+ m 2 H = 0. For m 2 < 0, we can define x — m 2 H l\m 2 \. Then the Eq. (|34l) can be written as 

2D,_ nrr ns _ ( e *Q 2 2 qV (l + l) {e*m 2 H f 1 (2 , *(1,^; 



a (w-0,T-0)- 2_, { _ 1 + {l+ 3 )x)2{ _ 1 + {l+ i )x)2 ~ h m 4^+ 2 X 3 X ) ( 58 ) 

The first singularity appears at x = 2 which is the superconductor-metal transition point and it yields a critical 
magnetic field H c = I ^r i . Expanding the conductivity around x = 2, we obtain for h — > 

1 h m- 1 fh 2 + h s 

(T X x [e*y qrj 



Pxx = — = 2*——(^^) (59) 



Hence, the leading term is h 2 which is consistent with the turn-on of the resistivity in the disorder-tuned transition. 
A comparison of this data and a best fit plot to the magnetic- field tuned data is shown in Fig. (JjJJ . The consilience 
of these results unifies the disorder and magnetic-field tuned transitions. Both have the same cause and their turn-on 
exponents resistivity as a function of the disorder and magnetic field are equivalent. Our predicition that v « 1/2 can 
certainly be tested further experimentally. In comparing with experiment, it is certainly best to know the width of 
the critical region. However, since heat capacity measurements are not available, we can only obtain a lower bound 
using the standard Ginzburg criterion. In our work, we calculated the conductivity in the zero temperature limin 
terms of the inverse correlation length, m 2 as, which to leading order in the field is given by 

~9, , 9 e*H , . e*H e*H c (H-H c ) e*H Cl , . 

m 2 (T - 0) = m 2 + — = a -T c + — = - cV - c> = -— -h 60 
nc nc he H c nc 

Substituting this expression into the Ginzburg criterion and assuming that the specific heat term is of order unity, we 
find that h = (H — H c )/H c w 10~ 6 sets the lower bound for the width of the critical region. 
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FIG. 5: The turn-on of the resistivity as a function of magnetic fiel d fo r MoGe near the superconductor-metal transition on 
the metallic side. The theory (slid line) curve is determined by Eq. (1591 . The best-fit curve, 70Q(H — H C ) 2A17 is in excellent 
agreement with the theory and with the experimental result of R C (H — -ff c ) M , 1 < fJ, < 3. The data were taken from N. Mason's 



III. GLOBAL SUMMARY 



All the bosonic models considered here are obtained from an action of the form, 

2l 1 



= / d d x dr^T^r) 



to - V -A(x,t) 

cn 



+ Jd d xJ drJ2^n( x > T ^*U x > T ^=T 2 =r + U J dTj2K( x > T )*>> T tf 

a a 

+77 J d d x I dr 1 dT 2 Y / ^,r 1 )\^-\^(x,r 2 )\ Tl=T2=T + (3 J d d x J dr^^^^^r)}* . (61) 

a ~ * * a, b 

The parameters which determine the conductivity are the dissipation to the Edwards- Anderson order parameter 
q ab , the quartic interaction U, temperature /3 = X/KbT and magnetic field B. In the following, we provide a brief 
table which catalogues all the distinct cases that arise and the dependence of the conductivity on the above factors. 
The table exhausts all known Bose models in which the conductivity has or can be calculated. 



IV. FINAL REMARKS 



We have shown here that the bosonic conductivity in the phase glass is robust to both interactions and a non-zero 
magnetic field. Further, we find complete agreement with between theory and experiment for the disorder-tuned and 
the magnetic-field tuned turn-on of the resistivity in 2D films. In general, the resistivity turns on as 

p{T = 0)<xp c (g~g c ) 2zv (62) 

where g is either disorder or magnetic field. The significance that a single expression with universal exponents 
describes both the disorder-tuned data as well as the field-tuned data is two-fold: 1) both the disorder and magnetic- 
field transitions share a common origin and 2) both are driven by superconducting fluctuations, that is, both phases are 
bosonic rather than fermionic in character. In fact, the power-law turn-on of the resistivity places severe restrictions 
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TABLE I: Conductivity for bosonic models as a function of system paraments. 



Distinct Cases 



7] = q = 



77/O q = 



-3- 



b = 
B = 



r = 

T = 



17 = 
17/0 



I 
I 



I 

I(QD) 00 SC(RC)ii 



M 



e* 4gy /-1 _i_ 
h 3^\ L 



,(7 
7rm J 



M 



B = 


T/0 


f/ = 


00 PC 18 


f£[^e-/ T + (fS) 3 ] (T«m) 

2^1n(^)(r;«m;m«T) 

^( ?7 »m;m«T^. 


Box left+Box above 6 


B = 


T ^ 


U ± 


00 so^a 


(7 renormalize m^ a 


U and q renormalize m a 


B ^ 


T = 


[7 = 


1 


I 




B / 


T = 


(7/0 


1 


T or no SC 




B/0 


T/0 


[7 = 


00 PC 


^[f^ /T + (^) 3 ] 

71/ 1 '/ OTTO ' J 

(m|f «mT);T« m) 

4e 2 r8r)Tm 2 -rh/T . / tt^T n3i 
hi? I- + \ 3ft 3 7 J 

(jjm < m|f « mT; T«m) 
2e2T RnC T 1 + llm *H 

h»7 WmJ ~ 32m 2 J 

(m < T;?7 < rh,m H ) 

e 2 r,T , m H] 
hm 2 1 ' 2m 2 1 

(m < 7/ < m« < T)^a 


Box left+Box above 6 


B / 


T^O 


[7^0 


00 sc 


L7 and B renormalize m a 


U,B and g renormalize m a 



a U,B and q will renormalize m and change the dependence of temperature T thus change the behavior of the conductivity indirectly. 
Sec scctiomRole of interaction. 
6 The conductivity in this case is equal to a (77 ^ 0, q = 0,T + 0,U = 0) plus a(r) ^ 0, 9 + 0,T = 0,C7 = 0). 

C PC, SC, I, M.QD and RC represent perfect conductor, superconductor, insulator, metal, quantum disorder regime and renormalized 
classical regime respectively. 



any permissible theory of the Bose metal. Further, identical results hold in D=2 as well as in D=3. Hence, the 
phase glass is a candidate to explain the experimental findings of the leveling of the resistivity as T — > in thin 
films that nominally display a superconductor-insulator transition. The fact that the metallic phase persists in D=3 
is particularly revealing in light of the recent experiments which find that in the putatative vortex glass phase, the 
I — V characteristics remain linear as opposed to the non-linear behaviour anticipated for a true superconducting 
phase. As it turns out, no explicit calculation of the conductivity has ever been performed for the vortex glass. The 
only difference between the vortex glass and the phase glass considered here is that in the former, the disorder in 
the Jij's is bounded in the interval [—1,1]. This is not a substantive difference. Hence, the model considered here 
contains the complexity of the vortex glass and hence should suffice to describe this model. Further, as we have 
shown above and elsewhereiS, coupling the electromagnetic gauge directly to the spin glass order parameter does not 
proves conclusively that the stiffness vanishes. In light of this result and the results derived here, we conclude that 
the vortex glass is in fact a metal. Ultimately, if our proposal is correct, the metallic phase in the thin films has been 
seen before. Namely, it is the vortex glass. 



V. APPENDIX A: 

The propagator, Eq. (J2J, describes bosonic transport for the free energy given by Eq. 10: 

^"gauss = {k 2 +ul + v\u n \+m 2 )\il; a (k,uj n )\ 2 

-f3q <5 <1>B ,oV°(k,w n )[V 6 (k,w n )]*, (63) 

a,6,k,u> n 

Consequently, in replica space, the inverse propagator is 

Gjjfcw,,) = (m a + k 2 + u 2 + V \lo\ -pq ab 5^Y a l (64) 

That is, 

J2 (k, w n ) (m 2 + k 2 + lo 2 + v \u\ - pq bc 6 Ut0 ) = l ac (65) 

6 
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To invert the matrix, we set 

G° b (k, w„) = A(k, Lj n )5 ab + B(k, u n )Pq ab (66) 

and 

C(k,w n ) = (m 2 +k 2 +w 2 +?7|w|) (67) 
and then substitute into Eq. (|64|l to obtain 

l ac = ^(^^6^ + B(k,uj n )(3q ab ) (C(k,uj n )S bc - (3q bc S ufi ) 

b 

= A(k, w„)C(k, ^„)5 ac + (B(k, cj„)C(k, «„) - A(k, w n )<5 w , ) /V c + B(k, u n )6 ufi E 9°V C - ( 68 ) 

b 

The last term is a product of Parisi matrices which is parameterized by the function q and q(x) (0 < x < 1), where 
q and (7(2;) are the diagonal and off-diagonal elements of q in the continous representation. The sum over the replica 
index b in the last term in Eq. (|68|l will be proportional to n and hence vanish in the n — > limit. Consequently, we 
simplify Eq. (|68|) to obtain 

A(k,w n )C(k,w n ) = 1 
B(k, w„)C(k, « n ) - A(k, u n )S u , = 0. (69) 

The solution to these equations is 

A(k,u n ) = — — ^ — - = 2TTTT 2 i — n = G (k,cj n ) 
C(k, u„) + k z + w z + 77|w| 

= c(k,o, n ) = RTW = Go(k ' w " ) ' 5 ^ - (70) 

Eq.© is then derived. 

VI. APPENDIX B: 



The point of this Appendix is to show that even if an Ohmic dissipation term is put in by hand as done recently 29 , 
nothing of a substantive nature changes. Hence, the claim that D cx l^l 1 / 2 leads to a superconducting state instead 
of a metallic one is false. To avoid any confusion, we redo the cumulant expansion. 

The effective action for our problem is given by 



+ yEE/ f dTdT '\ E ^{i^-ofl.Vj-CW-a^C^-^W-a^tr'))]} 

a,b (ij) J ° J ° a=+l,-l 

+C.C. 

dr ex p K Q - ° h M) - A v) + c - c - 



Jo 
2 



(71) 



a i] 



Using the Hubbard-Stratanovich transformation and introducing the order parameters ^"(fc,r) and Q a ^ v (k, t\, T2), 
we have the effective free energy^, 



F eS [V,Q] = ^2[ dr*«(fc,r)[*J(fe ) r)]*+ E / dr / ^' Q^(fc,T,r')[Q^(^r,T')] a 

I JO L U Uf JO JO 

" a.b.k,k' 

f exp Uj^drJ2 y/Jo~(k)^(k, r)S°(-k, t) 



16 



V f dr f dT'^2J{k)J{k')Qf u {k 1 k', r, r')S;(k, r)S b v (k', r') + F dls ) ) 

a,b,k,k' J ° J ° J/ Q 



dis 



dr / dr'- 



7r(r — t' 



1- ]T S;(k,r)S^(-k,T') 

a,fj,,k 



where J(k) = J^2 R cos(kR) 1 Jo(k) — JoJ2 R cos(kR) (R is the displacement vector of nearest neighbour) and 



(72) 



A) 



(73) 



Here, Zq = Tr[exp(— /3Hq)]. The Hamiltonian Hq and Hq are denned as, 



Hn = - 



—r(- 



H, 



eff _ 



t-y(- 



(74) 



where a is the index of replica and i is the site label. 

The only difference from ReS^ with respect to the cumulant expansion is the dissipation term. However, inspection 
of the free energy indicates that this term can be combined with the second-to-last term by defining a new Q matrix, 



V2J(k)J(k')Q^(k : k',r,T')^S ab S, u 6 k+k , fi . V m y/2J(k)J(k')Qf v (k,k',T,T') 



(75) 



The new free energy, 



f exp ^ J* drr^yfiW*%{Jz,T)Sl(-k,T) 



E 



, ^0 



dr / dr' 



^2J{k)J{k')Qf v {k, k', r, r') S a Jk, r)S b u (k', r') 



(76) 



is identical with Q replaced with Q. Hence, the solution to Q is the same as that for Q. Once the transformation, Eq. 
(17511 is undone, we recover that the Ohmic dissipation term simply renormalizcs the total coefficient of \uj\, leaving 
the glass and metal intact, contrary to the claims of Ikeda2^. 
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